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Abstract 

Almost all normal-incommensurate phase transitions observed experimen- 
tally are continuous. We show that there is not any theoretical reason for 
this general behaviour in perfect crystals. A normal-incommensurate phase 
transition that is not too far from the mean-field tricritical point should 
be discontinuous and it is highly improbable that so far reported normal- 
incommensurate phase transitions lie very far from this point. To under- 
stand this behaviour we study influence of defects on a hypothetical first-order 
normal-incommensurate phase transition in a pure material. We have found 
that this influence is strikingly different from that on other kinds of first-order 
phase transitions. The change of the discontinuity of the order parameter at 
the transition is negative and formally diverges within our approximate the- 
ory. At the same time the diminishing of the phase transition temperature 
remains finite. We interpret these results as an indication that at least some 
of the observed seemingly second-order normal-incommensurate transitions 

would be first-order transitions in defectless crystals. 
PACS numbers: 61.44.Fw, 64.60.-i 
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I. INTRODUCTION 



Almost all observed structural normal-incommensurate (N-IC) phase transitions are con- 
tinuous (see for instance Ref. [I] as example of continuous cases and Ref. [|for rare examples 
of discontinuous ones). This is in contrast with other types of structural phase transitions 
that are generally discontinuous. Such a discontinuity is not surprising because, as was 
pointed out by many authors!, the second-order (continuous) phase transitions convert into 
first-order (discontinuous) ones because of many factors like shear deformations, crystal 
anisotropy, etc., except for phase transitions involving uniaxial ferroelectrics and ferroelas- 
tics. These factors are operative by virtue of order parameter fluctuations. We show that 
also a broad class of normal-incommensurate (N-IC) phase transitions should be discon- 
tinuous in pure materials due to the solid-state elasticity. We refer to transitions that are 
not too far from the mean- field tri critical point. Then a question arises: why all observed 
experimentally N-IC phase transition are continuous? 

In perfect crystals the N-IC first-order phase transitions should exist but it is reasonable 
to assume that in real crystals the defects can smear them@. Indeed, it is well known! that 
some defects influence dramatically on the incommensurate phases, destroying the long- 
range order at any concentration. Therefore strictly speaking, there are no N-IC phase 
transitions in real crystals or in other words, these transitions are always smeared. However, 
this smearing is not observed in the experiments and N-IC phase transitions appears usually 
as "ordinary" transitions, i.e. without apparent smearing in the measured anomalies. As we 
have mentioned, what is surprising is that almost all N-IC phase transitions are second-order 
ones. We try to understand this general behaviour showing that first-order N-IC transitions 
can assume a continuous appearance if the crystal contains defects. 

The paper is organized as follows. In Sec. II we demonstrate that close to the mean-field 
tricritical point a second-order N-IC phase transition converts into a first-order one because 
of the combined effect of the solid-state elasticity and the order parameter fluctuations. 
It is important that the conversion takes place even in the framework of the first order 
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perturbation theory of effects of critical fluctuations. The vicinity of the tricritical point 
where this conversion takes place is shown to be broad. In Sec. Ill we present the method to 
consider influence of defects on first-order phase transitionsiB within first order perturbation 
theory where only linear terms in defect concentration taken into account (see e.g. Ref. §,0). 
To simplify the treatment, the phase transition in the pure material is supposed to be 
describable within the classical Landau theory. In Sec. IV we apply the method to a first- 
order N-IC phase transition. Because of the specific features of a incommensurate phase, 
this theory is not rigorously applicable here but we use the fact that N-IC transitions can 
be considered as a limiting case of some normal- commensurate (N-C) transitions. So we 
consider influence of defects on a normal-commensurate (N-C) first-order phase transition 
where the perturbation theory is valid, and then trace the changes produced by defects when 
the phase transition approaches a N-IC one. As we will see, we cannot reach the N-IC phase 
transition because that means to go beyond the applicability of the perturbation theory, but 
we can reveal some qualitative trends which provide us a base for the speculations. In Sec.V 
we show that the same trends take place when the first-order phase transition in the pure 
material is due to the combined effect of order parameter fluctuations and the solid-state 
elasticity as in Sec. II. Finally, the Sec. VI is devoted to the conclusions. 



II. FLUCTUATION-INDUCED FIRST-ORDER N-IC PHASE TRANSITIONS 



Among several factors that are relevant in the conversion of a second-order phase tran- 
sition into a first-order one the most important one seems to be the solid-state elasticity. 



Indeed, as it is shown in Ref. [10], close to the mean-field tricritical point the combined effect 
of the order parameter fluctuations and the solid-state elasticity is so strong that the second- 
order phase transition converts into a first-order one quite far away from the region where 
the critical fluctuations become too much important. In other words, this converted phase 
transition takes place in the region of applicability of the Landau theory where the effects of 
the critical fluctuations can be considered within the first order perturbation theory. This 



question was reconsidered for structural displacive transitions in Ref. [11], and it was shown 
that the region where the above-mentioned mechanism of conversion of a second-order into 
first-order phase transition results operative is very broad. Only the case of one component 
order parameter was considered. We shall show now how the same conclusions are valid for 
a N-IC phase transition. 

We restrict ourselves to the case of one-k incommensurate phase, i.e. the case of two 
component order parameter (designated as rj(r) = (771 (r), 7y 2 (r)))& Then, the Landau 
continuous-media potential takes the form: 

u ik ) = J [ip v {v) + Vu(rj, u ik )} dv (1) 



with 



<Pv(ri) = ^i(r) + % 2 (r)) + Jb(t£ (r) + % 2 (r)) 2 + 

+ lc(vl(r) + % 2 (r)) 3 + ^((Vr^r)) 2 + (Vr^r)) 2 ) (2) 



and 



X 1 

<p u (ri,u) = r(r/ 2 (r) + r/ 2 (r))u /z (r) + -Ku 2 u (r) + n{u ik (r) - -u u (r)5 ik ) 2 (3) 

where K is the bulk modulus and /1 is the shear modulus (to compare with Ref. 0). We 
shall assume that the coefficient A in Eq. ^| is the only temperature dependent coefficient: 
A = a(T — T c ) where a is a positive constant and T c would be the second-order transi- 
tion temperature if coupling between elasticity and fluctuations is neglected. The other 
coefficients may depend on pressure but we consider them positive constants. 

As usually, it is assumed that the Landau potential is obtained as result of a partial 
integration of the partition function over all degrees of freedom, except those that correspond 
to the long- wave Fourier components of 77 (r) and tt«;(r). To find thermodynamic quantities, 
one has to integrate besides all the degrees of freedom except those probed in experiment. 
For our purpose it is enough to leave unintegrated the zero Fourier component of the order 



parameter. As in Ref. |TT] we minimize over the elastic degrees of freedom. At this point one 



has to make out the spatially homogeneous (uf^) and inhomogeneous (ufj) parts of strain: 

uy(r) = + u§ = «g> + W [WW + k jU \k)} e* kr (4) 

where u(k) is k-Fourier component of the displacement vector. Minimizing over the homo- 
geneous and the inhomogeneous strain separately one obtains: 
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J if u (r], U)dv = -7tf((^L falk^l-k + ^2k^2-k]^) - 



r 2 

J2 [»7iki»7i-ki-k + »72ki»72-ki-k] [^i- k2 ^ik 2 +k + ^2-k 2 ^2k 2 +k] (5) 

where 7 = X + |/x and we put the system volume equal to unity. 

Taking these terms into account, the Landau potential can be rewritten (limiting us to 
harmonic terms) as: 

$(77 ,.., 77 k ,...) = $o(T7 ) + ^ h {rj (6) 

where $0 is the Landau potential that depends only on the zero Fourier components of rj 
(designated for simplicity as rj = (771,772)): 

MVo) = \A{vl + ril) + \m + rjtf + ±C(tf + r£? (7) 
and $^ are the harmonic terms of the Landau potential that depend on r] and r/ k _^ : 

^h(»7o.»7k^o) = SiUkM^-k (8) 

k^0 Z 

where 

Mn = A + 35i77 2 + 5 % 2 + C(5rtf + 77* + 6t7 2 77 2 ) + -DA; 2 (9) 
M 22 = A + ZBwl + S77 2 + C(5^ + 77^ + 6t7 2 77 2 ) + Dk 2 (10) 
Mi 2 = (2£ + 3A)77!77 2 + 4C(t7i77 2 3 + Tjfo) (11) 



In the preceding equations the constants B, B\ are defined as: 

2r 2 



(12) 



_ 2r 2 4r 2 ~ 16r 2 /z 

Bi = B — = B H 

37 97^7 



5 + A 



(13) 



As we have mentioned, we must integrate over the degrees of freedom corresponding to 
the non-zero Fourier components of the order parameter according to the Gibbs distribution. 
Thus, the incomplete thermodynamic potential is: 



$ = $o(r, ) - Tin {/ exp n^d** 
As is quadratic in 771^ and 772^, the integration can be easily performed, obtaining: 

$(T7 ) = $oM + \ |E In (M U M 22 - M 2 2 ) | 



(14) 



(15) 



Using the fact that in the incommensurate phase, the system is invariant respect to the 
phase of the order parameter vector 77, we can choose it freely. Taking 1)1 ^ 0, 1)2 = we 
considerably simplify the expression of the thermodynamic potential because in this case 
M22 = 0. Then, carrying out the integration of the second term of Eq. [15], and bearing in 
mind that Dk^ > A + SB^j + 5Cr]f , we obtain: 



*(»7o) = g 



A + 



k T 

""max- 1 

2n 2 D 



(35i + B) 



B + 



k T 

max - 1 - 



ll,i,x *-6C 



7T 2 D 



T 



r2nm i{(A + 3B^ + 5C V ^ 



(A + Bri + Crt? 12 ) 



(16) 



There is a trivial renormalization of the reference energy taken into account and denoted 
by the symbol $ in the left side of the Eq. |TB[ From the first two terms, one can see 
that the fluctuations provoke a renormalization also in the coefficients A and B. This 
renormalization also takes place in the rest of the terms because of higher order corrections. 
Since this renormalization is not of interest within our phenomenological theory, henceforth 
that all the coefficients correspond to their renormalizated values. Then, the final form 



of the thermodynamic potential including the effects of order parameter fluctuations and 
solid-state elasticity is: 

I'M = l -Arf x + \Br(i + ^Cvf - 



-m 



A + 3B lV j + 5C V tf 2 + (A + B V l + Crtf' 2 } (17) 



where m = T/12ttD 3/2 . 



The same discussion as in Ref. |TT] about the final character of the phase transition can 
be carried out now for a N-IC one. As in that work we can see that close to the mean-field 
tricritical point, the order parameter at the transition temperature is such that: 

B lV 2 ^A, Crjf (18) 

As in Ref. [11] we find the value of the coefficient A at the transition temperature (A') and 



the order parameter discontinuity at the transition (f] f2 ): 

" ? = (1f) 2/3(3Bi) (20) 

We must conclude that the combined effect of order parameter fluctuations and the solid- 
state elasticity over a N-IC phase transition forces it to be discontinuous. This assertion 
is correct in a region close to the mean-field tricritical point, i.e. as can be see in Ref. |ll] 
in a region such that B < fiB/K. It is highly unlikely that all the reported N-IC phase 
transitions lie out of this region because it is very broad for structural systems. Such a 
conclusion is only applicable when dealing with perfect crystals because in our treatment 
we have forgotten about the influence of defects. 



III. INFLUENCE OF DEFECTS: THE METHOD 



In this section we expose briefly the simplest version of treatment of influence of defects on 
first-order phase transitions^. This allows us to present formulas that are easy to generalize 
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in the cases of our interest. What is more important: in these formulas we keep a term that, 
although omitted in previous papers, it is essential in the context of this work. 

Within this section we consider one component order parameter (rf) and "random local 
field" (RLF) defectsllll!. The Landau continuous-media thermodynamic potential without 
defects has the form: 

$ = J (p(rj)dv (21) 

with 

<p{rj) = \Arf{v) + iiV(r) + l -C V \v) + ^(Vr?(r)) 2 (22) 

where A = a(T — T s ), T s is the temperature at the spinode of the symmetrical phase and 
all other coefficients are supposed to be temperature- independent. As we are considering a 
first-order transition we assume that B < and all the other coefficients positive. 

The first-order transition temperature and the value of the order parameter discontinuity 
at the transition is obtained with V?y(r) = and using the equilibrium equations: dQ/drj = 
(minimum of thermodynamic potential condition) and $ = (equality of thermodynamic 
potentials of the two phases condition). One has: 

A + B?) 2 + C?j 4 = (23) 



Arj 2 + \b^ + \cf = (24) 

where A and rf 2 are the values corresponding to the first-order phase transition temperature. 
From these equations: 

f - ~ (26) 

in a perfect crystal. 



To take into account RLF defects one has to add to the thermodynamic potential density 
<p(rj) in Eq. |22] a term: 

- h{r)r}{r) (27) 
where h(r) is a "random local field" in the form: 

Mr)=]TV(r-r,) (28) 

3 

where hj = ±if, if is a coefficient characterizing the "strength" of the defect localized at 
r = rj, and the sign of hj is random. 

The order parameter decomposition into Fourier (r/k) components is the basic tool to 
study the influence of defects, so expressing: 

77(r)=5> k e ikr (29) 

k 

we consider the thermodynamic potential $ as a function of all Fourier components. The 
next step is to minimize this function over all Fourier components except for the zero com- 
ponent (k = 0, designated by 770) that at this stage is consider as a fixed parameter. One 
obtains to the linear approximation: 

= A + Wrft + hC4 + Dk 2 (30) 

Then, we use these equilibrium values in the thermodynamic potential and average over the 
defect positions obtaining! I: 

Hvo) = \Avl + \Bvt + + ^j-M + Wrjl + 5C V ^ (31) 



where N is the defect concentration. Instead of Eqs. ^3] and |24] that give to us the values of 
the order parameter discontinuity at the transition and the transition temperature, one has 
now: 

(32) 



AT) 2 + \bt& + \c?f + 2n{A + 3Brf + 5C^) 1/2 = 2nA^ 2 (33) 

where n = NH 2 /8ttD 3 ^ 2 , A and rj 2 are the values at the first-order phase transition temper- 
ature. Let us mention that in previous papersBiB the right hand side of Eq. |33| correspond- 
ing to thermodynamic potential of the symmetrical phase was set equal to zero. This term 
would not really modify the previous papers results but plays a basic role for a N-IC phase 
transition. 

The applicability conditions of these formulas at first-order phase transition temperature 
can be obtained in the same way as in Ref. and implies for the symmetric and non- 
symmetric phases: 

C 3/2 

n^r < 1 ( 34 ) 

where the numerical factor (of magnitude order of unity) has been omitted. 

Designating the changes produced by defects in A and rj 2 (Eqs. [25| and ^6|) as 5 A and 
Srj 2 , one can linearize Eqs. [32] and |33| with respect to this values to obtain: 

8A = -2nM (35) 



-\An\ 



l_C_ 

w 2 



(36) 



One sees that RLF defects provoke both lowering of the phase transition temperature and 
diminishing of the discontinuity of the order parameter at the transition^ According to 
Eq. we cannot set B = in the last equation. Moreover this linear approximation to 



Eqs. [32] and [33] lose their applicability at the same time that the equations do. So nonlinear 
corrections to equations |35| and |36| would be inconsistent with the applicability criterion for 



equations |32j and 
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IV. APPROACHING A NORMAL-INCOMMENSURATE TRANSITION: 
DEFECTS AND LANDAU THEORY 



In as follows we consider N-IC phase transition as a special case of phase transition with 
two component order parameter (rji and 772) corresponding to isotropy in order parameter 
space. This means that all invariants that we must consider in Landau thermodynamics 
potential are power of the second order invariant^. In other words, taking the Landau 
potential density in the form : 

ip{ri) = l -A{r,l{v) + r, 2 2 (v)) + ifl^r) + 77 2 (r)) 2 + -^(r^r) + 

+~C(7? 2 (r) + r,l{v)f + ~D((V V i(r)T + (Vr? 2 (r)) 2 ) (37) 

we approach a N-IC first-order transition if i?2 — ^ and all other coefficients are considered, 
positive except Bi and temperature independent except A = a(T — T s ). Now to make 
allowances for defects we add to the thermodynamic potential density a term: 

- /ii(r)7?i(r) - h 2 (r)r] 2 (r) (38) 

where h\ and h 2 are the two components of a "random local field": 

hi, 2 (v) = J2h{Ar-T j ) (39) 

j 

We assume that ((h{) 2 } = ((hi) 2 } = H 2 and that the "direction" of the "vector" h = (hi, h 2 ) 
is random, i.e. (h{h 2 ) = 0, where ( ) designates the average. 

We apply the method of the preceding section to find out the thermodynamic potential as 
a function of the zero Fourier components (k = 0, designated for simplicity by 771 and 772) . It 
is convenient to suppose that B 2 > 0. The same results can be obtained for the case B 2 < 
but the expressions are considerably longer in the latter case. Thus in what follows we take 
the case in what expressions are straightforward, i.e. B 2 > 0. Then, non symmetrical phase 
equilibrium state corresponds to rji 7^ 0, rj 2 = or vice versa. Taking rj 2 = we obtain 
instead of equation [TTJ for a one component order parameter, the thermodynamic potential 
in the form: 
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r) 2 ) = X -A^\ + -B lV f + -Cr,l + 

+n {(A + 35^ + 5C7/?) 1 / 2 + (A + + B^ 2 + CV^) 1 / 2 } (40) 

and instead of equilibrium equations |23] and |24] : 

(41) 

M)\ + -B^ + icfTf + 2n {(A + 3S^ 2 + 5C^) 1 / 2 + 

+ (I + (B 1 + £? 2 )f/ 2 + Cr,t) l/2 } = AnA 1 ' 2 (42) 

where as before n = NH 2 /8ttD 3 ^ 2 , A and rj 2 are the values corresponding to the first-order 
phase transition temperature. 

The applicability of these equations at the first-order phase transition for the non- 
symmetric phase in the limit B 2 —>■ (to try to describe the N-IC phase transition), implies: 

£3/2 

%l/2 R 3/2 < 1 ( 43 ) 

revealing only the dependence on thermodynamic potential coefficients. It is obvious that 
this condition is violated if B 2 is sufficiently small, and therefore we can approach but not 
reach the N-IC phase transition. 

In the same way as before the changes produced by defects in A and rff (Eqs. ^ and ^6|) . 
now for a two component order parameter potential, are given in the linear approximation 
by: 



B 2 C 



SA—AnJ^ (44) 



^-HM 9 ~™ + ^] (45) 



The applicability criterion of this linear approximation to Eqs. f|T] and [42| is found easily 
imposing the natural condition Srj 2 < rj 2 . In the limit B 2 —>■ 0: 



li= -Hmr\ (46) 
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= - 2n mWYK (47) 

and without the numerical factor, his applicability criterion is: 

^3/2 

%3/2 R l/2 < 1 ( 48 ) 
-°1 #2 

In this case, this applicability criterion is less restrictive than Eq. ^3|. This means that 
one can, in principle, go beyond the linear approximation remaining within the region of 
applicability of Eqs. f|l] and ^2] but it is not necessary for our discussion. 



As we can see from equations [46] and |47|, for an approach a N-IC transition, the correction 
to the transition temperature goes to zero while the correction to the order parameter 
discontinuity is negative and diverges. This means, of course, that the perturbation theory 
is no more applicable here. But also indicate that while defects suppress very strongly 
the order parameter discontinuity they do not influence (or weakly influence) the phase 
transition temperature. 



V. APPROACHING A NORMAL-INCOMMENSURATE TRANSITION: 
DEFECTS AND CRITICAL FLUCTUATIONS 

In the preceding section we have considered phase transitions that are first-order irrespec- 
tively to the order parameter fluctuations. However there exist another type of first-order 
phase transitions: those that are discontinuous due to critical fluctuations. In Sec. II we 
shown that structural N-IC phase transitions belongs to this class if it takes place in a broad 
vicinity of the mean- field tricritical point. In this section we study influence of defects on 
this type of transitions. 

We study a N-C phase transition of the same symmetry class that in Sec. IV and we 
follow the evolution of the results when a N-IC phase transition is approached. For the 
thermodynamic potential one has now: 
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$(77 , 77,, ...) = l -Aril + ±B lV t + l -C V t 



-m 



A + 3(3 V 2 + ZCtff 2 + (A+ (B 1 + 5 2 )^ + C^f' 2 } + 



^l, *0»7ik»7i-k + X^X 2 H^i) k)r] 2k r] 2 _ k (49) 

k k 



where as before m = T/12tiD 3 ^ 2 , but now making the definitions: 



~ 2r 2 

S x = S X - — (50) 



O r 2 4^,2 

^ B '-37f-3T" Bl + A <51) 



The coefficients 



X^fai, k) = A + 3/fy 2 + 5Ctf + Dk 2 - 3m {(A + 3(3^ + 5C^)V 2 (3/3 + 10Cr/ 2 ) + 

+ (A + (Si + S 2 )77? + Crt) 1 ' 2 ^ + B 2 + 2Cr,l)} (52) 

Xj 1 ^!, fc) = A + (Bi + B 2 )rji + Cfy* + DA; 2 - 3m {(A + 3/?r/ 2 + 5C^) 1/2 (3/3 + 10Cfy?)+ 
+ (A + (B, + B 2 )t£ + Ciif' 2 {B x + 5 2 + 2CV/ 2 )} (53) 

are calculated as in Sec. II performing the integration in Eq. [TJI, but now keeping out of this 
integration r] and rj^, where k' is a fixed wave vector. 

Defects are considered adding to this potential the term: 

- ^2(hiU]i- k + /i 2 k^2-k) (54) 

k 

Proceeding as before with this thermodynamic potential we obtain in the limit B 2 —>■ : 

/ 2(7 \ 1/3 

SA> = -2n(-) (55) 



< = -n ) [ ^ ) (56) 



2Ty :! /^ x ' 2 
V3mJ V3/?S 2 

close to the mean- field tricritical point. We see that for an approach a N-IC transition, the 
correction to the order parameter discontinuity is negative and diverges in the same way that 
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it do for the Landau theory (Sec. IV). Here the correction to the transition temperature also 
remains finite but do not tends to zero. As well as in the preceding section this indicates 
that while the defects suppress very strongly the order parameter discontinuity even for 
small defect concentrations but they weakly influence the phase transition temperature. 

VI. CONCLUSIONS 

Hypothetical N-IC first-order phase transitions are strikingly different from other "or- 
dinary" first-order phase transitions where phase temperature always diminishes because 
of influence of defects, see Fig.[l|. Lowering of transition temperature in Fig.|l| (a), (b) is 
consistent with the expectation that at some high enough defects concentration the phase 
transition temperature becomes equal to zero and the phase transition no more exists. The 
practical absence of this shift for a N-IC first-order phase transition (Fig .[I] (c)) seems to 
indicate that the crystal with defects "remembers" the first-order transition temperature 
in the pure material. This temperature might correspond to fairly rapid but still continu- 
ous changes in the satellite intensities with the temperature that might be interpreted as 
second-order transition. The results of this paper emphasize once more the strong influence 
of defects in the properties of incommensurate phases. 

We would like to thank S.A. Minyukov, J.J. Saenz and the members of his group specially 
to A. Garcia-Martin and L.S. Froufe for useful discussions. Special thanks to D. Sanchez- 
Soria for her help and support. 
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FIGURES 

FIG. 1. (a)Influence of "random temperature" defects on the phase diagram of a first-order 
phase transition with elasticity taken account. (b)Influence of "RLF" defects on the phase diagram 
of a first-order phase transition with elasticity taken account. If the elasticity is neglected the same 
behaviour that in (a) is obtained. (c)Expected influence of "RLF" defects on the phase diagram 
of a N-IC first-order phase transition 
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